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ABSTRACT 

Recent observations point to the presence of structured dust grains in the discs surround- 
ing young brown dwarfs, thus implying that the first stages of planet formation take place 
also in the sub-stellar regime. Here, we investigate the potential for planet formation around 
brown dwarfs and very low mass stars according to the sequential core accretion model of 
planet formation. We find that, for a brown dwarfs mass 0.05 M©, our models predict a max- 
imum planetary mass of ~ 5M e , orbiting with semi-major axis ~ 1AU. However, we note 
that the predictions for the mass - semi-major axis distribution are strongly dependent upon 
the models chosen for the disc surface density profiles and the assumed distribution of disc 
masses. In particular, if brown dwarf disc masses are of the order of a few Jupiter masses, 
Earth-mass planets might be relatively frequent, while if typical disc masses are only a frac- 
tion of Jupiter mass, we predict that planet formation would be extremely rare in the substellar 
regime. As the observational constraints on disc profiles, mass dependencies and their distri- 
butions are poor in the brown dwarf regime, we advise caution in validating theoretical models 
only on stars similar to the Sun and emphasise the need for observational data on planetary 
systems around a wide range of stellar masses. We also find that, unlike the situation around 
solar-like stars, Type-II migration is totally absent from the planet formation process around 
brown dwarfs, suggesting that any future observations of planets around brown dwarfs would 
provide a direct measure of the role of other types of migration. 
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1 INTRODUCTION 

Observations of planets proceed apace: currently there are 232 ex- 
tra solar planets known in 199 different planetary systemfl As is 
well known, many, if no t all, of these s ystems are radically differ- 
ent to the Solar System dBeer et al]| 20041) . often posing challenges 
to theories of planet formation. As time has progressed, a num- 
ber of planets have begun to be observed around low mass stars 
(M* ~ 0.2 - 0.3M o ), thus begging the natural question of whether 
this trend will continue and whether similar planetary systems are 
to be expected at even lower masses, down in the brown dwarf 
regime. 

Observations suggest that there does not seem to be any dis- 
continuity between the properties of the discs arou nd stars and 
those around brown dwarfs: Spitzer observations by Luhm anet al.l 
j2005h found that disc fractions for stars and brown dwarfs in the IC 
348 and Chamaeleon I cl usters were similar, thus suggesting a sim- 
ilar disc lifet i me (se e also lSterzik et al. 2004). The measurement by 
IScholz et alj d2006h of the Taurus star-forming region indicates that 
many of the discs around brown dwarfs must be relatively large in 
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size, extending beyond 10AU, and have masses up to 12% of the 
central brown dwarf, with the average being 2.6%, compared with 
an average value of 1.9% for low mass main-sequence stars (note, 
however, that such disc mass estimates are very uncertain). Hence 
it appears that discs around brown dwarfs have similar properties 
and are as ubiquitous as discs around main sequence stars. 

Observations indicate that grain growth, the precursor to 
planet formation in the sequential core accretion model , can occur 
efficie ntly in the small discs surrounding brown dwarfs dApai et all 
2005). This then poses such questions as: what is the potential 
to form planets around very low mass stars; what would be their 
mass-semi-major axis distribution and would they be detectable? If 
formed or forming planets could be observed around brown dwarfs, 
these would provide highly valuable, low mass calibration points 
for the assessment of planet formation models over a wide range of 
stellar masses. In view of this, and in view of forthcoming missions 
such as the Terrestrial Planet Finder and Darwin, aimed at look- 
ing for Earth like planets, it is important to estimate the likelihood 
of the presence of planets around stars of different masses, and in 
particular around brown dwarfs. 

Planet formation models currently fall into two categories: (i) 
the sequential core accretion model, which can account for the for- 
mation of both rocky planets and gas giants and (ii) the gravita- 
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tional instability model, which can only account for the formation 
of massive gas giants. Th e sequential core accretion model of planet 
formation I Polla ck et al ] fl996h assumes that planets form through 
a process in which first dust grains initially present in the pro- 
toplanetary disc condense to form kilometre-scale planetesimals; 
the collisional accumulation of planetesimals forms the basis for 
the runaway growth of a protop lanetary core of rock and/or ice 
jSafronovH l972; Wet herilll l 1980), and finally if a critical mass has 
been reached, the pr otoplanetary cores can begin to rapidly accrete 
a gaseous en velope dMizund[l98(il: iBodenheimer & PoUackll 19861 ; 
|PollacketalJll996h . The gravitational instability model does not 
need to form the initial rocky core, and assumes that the proto- 
planetary disc may become gravitationally unstable and hence form 
gas- giants via gravitational collapse of high density regions in the 
disc lBossll 998. 2000: iMaver et alj2004l) . Whilst it is possible that 
some of the observed planets formed via a gravitational-instability 
scenario (and indeed also planetary mass compan ions to brown 
dwarfs. See, for example, the case of 2M AS S 1 207B . | Chauvin et all 
120051 : lLodato et afll2005h . the results of iMatsuo et alj j2007l) sug- 
gest that around 90% of the observed extra-solar planets have a 
mass - semi-major axis relation consistent with the core-accretion 
model and that the metallicity dependence of the distribution of 
these planets can only be explained using the core-accretion model. 

In a seri es of papers, llda & Linl J2004al lbl |2005|, henceforth 
lD-.lt llL2k llL3b attempted to construct a semi-analytic model of the 
combined core accretion and migration process which would allow 
quantitative comparison to be made between the observed distri- 
bution of extra-solar planets and that predicted by the theoretical 
core accretion model. They concentrated on investigating the ob- 
served mass - semi-major axis distribution of extrasolar planets, 
shown here in Fig. Q] This model starts with kilometre-size plan- 
etesimals, embedded in a two component gas and dust disc. The 
subsequent core accretion of planetesimals onto the protoplanet is 
modeled using the combined results of analytic models and N-body 
simulations. Once the cores reach a given critical mass, gas accre- 
tion is assumed to start, this process being modeled using numerical 
results for Kelvin-Helmholtz contraction of the gas envelope. This 
model also uses a number of prescriptions to enforce the limita- 
tion to core growth caused by various effects, such as isolation and 
migration. 

Whilst this kind of modeling is admittedly oversimplified, the 
models of IL1-IL3 do succeed in reproducing several features of 
the observed distribution of extra-solar planets, such as the metal- 
licity dependence for the probability of planet formation, the de- 
pendancy on stellar mass of the distribution of planetary masses as 
well as giving realistic mass ranges for the outcome of the planet 
formation process. However, the models experience problems in 
reproducing observations of intermediate mass planets with semi- 
major axes less than 3 AU. It is possible that these difficulties stem 
from the neglect of Type I migration in these models - see sections 
|2,4| and [3~!2l for further consideration. 

Clearly a self-consistent disc processing and core growth 
mod el for multiple interacting planets would be more appropriate 
(e.g.lLaughlin et alj206H;|Papaloizou & Nelsonl2005l;lAlibert et al.l 
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But given that the models of IL1-IL3 are currently the only extant 
attempt at providing a general framework for examining the out- 
come of the planetary formation and migration mechanism, they 
provide a useful framework for the exploration of the extensive pa- 
rameter space and they allow us to start examining in a statistical 
manner the potential for planet formation in discs around brown 
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Figure 1. Mass - Semi-major axis distribution of all currently known extra- 
solar planets. 



dwarfs and the potential for practical observations in the near fu- 
ture. 

To address the issue of planet formation around brown dwarfs 
under the core accretion model, we therefore extend the core- 
accretion model of IL1-IL3 down to the low mass regime. Once the 
basic structure of the model is in place, it can be used to investigate 
the effect on the planetary mass - semi-major axis distribution of 
varying key physical parameters. These include: the surface den- 
sity profile of the disc, the size and mass of the disc, the rate of 
disc clearance, the effect of stellar mass and luminosity as well as 
migration. 



2 PLANET FORMATION MODEL 

The model we adopt here is pri mar il y bas ed on the model of se- 
quential accretion presented in IlL ll - \MJ& In this section we ini- 
tially briefly recount the principle components which pertain to our 
investigation before going on to discuss in more detail the refine- 
ments and additions which are required to extend the model down 
to the brown dwarf regime. Readers requi ring more detail on the 
original model construction should refer to IlL ll - llOl 



2.1 Disc Properties and Initial Conditions 

The protoplanetary disc is modeled using a two component surface- 
density profile, £ = S rf + Z s , to account for the dust and gas, with: 



(1) 
(2) 



/ a \ k (M*\' , 
^ = Wo^^j^jgcm- 



where M+ is the mass of the central star, a is the radial distance 
from it, and r\ ice is a step-function accounting for the increase in 
dust density at the snow line such that rj ice (a < a ice ) = 1 and r\ ice {a > 
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a ke) = 4.2. The radial location of the snow-li ne a ice is set by the 
stellar luminosity, L*, such that dHavashi|[l98ll) : 



lar discs on T Tauri stars dHartmann et al.ll998l:lcaivet et alj200oh 



AU. 



(3) 



The factors (M*/M )" allow an investigation of the depen- 
dence of the disc mass on stellar mass (here we examine the cases 
n = 1 and 2). The dissipation of the gaseous component is mod- 
eled as a simple exponential decay with characteristic timescale, 
Tdi K , chosen to be ~ 10 6 — 10 7 yrs, as sugge sted by observation s 
of the infrared excess in young stellar objects dHaisch et alJl200lh . 
The dust component is assumed to be time-independent. Setting 
the parameters f = 1 an d k = -3/2 wou ld give the minimum mass 
nebula model (MMNM) dHavashi|[l98lh . We investigate the effect 
of varying the a b solu te surface density of the disc by varying fo. 

While IlLlI - llL3l only consider the case where the surface den- 
sity profile is a power law with index k = -3/2, here we also inves- 
tigate the effects of having k = — 1 (see section [3~.4. ll for motivation 
and discussion). 

It should also be noted that the simple exponential decay of 
the gaseous disc is the only evolution of the disc to be considered. 
The solid component is taken to be constant and neither component 
is modeled self-consistently with the growing mass of the planet. 
A self-consistent disc processing and core growth m odel would 
clearly be more appropriate (e.g. lAlibert et al.ll2005l) . but would 
also be more computationally demanding, so we adopt the simpli- 
fied approach in this investigation. 

To investigate the formation of planets around brown dwarfs, 
it is clear that we need to extend the semi-analytic model down to 
M+ < 0.084M o . However, the model depends on M t in two main 
ways: (i) Direct dependence on M t via £ oc MJ; and (ii) Indirect 
dependence on M*, e.g. the dependence of a ice and hence on L*. 
The detailed effects of varying L+ are discussed below in sections 
l23land l3A3l 



2.2 Core Growth and Gas Accretion 

We follow the prescription for core growth and gas accretion as 
presented i nllLll with a core a ccretion timescale for a planet of 



mass M of dKokubo & Ida1l2002h 
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where m = 10 18 g is the assumed mass of the planetesimals accreted 
by t he planetary core. 

Irkomaetalj (2000) find that once a critical core mass of 



M cril = 10 



M 



v lO-^^yr- 1 , 

has been reached, gas accretion can commence at a rate of 



(5) 



(6) 



Once a planet has grown to such as size that it opens a gap in 
the disc, the gas accretion rate is regulated by the rate at which the 
disc can process material and supply it across the gap to the planet, 
use a mass transfer rate given by the inferred rate for protostel- 
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yr 
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However, once the gap is opened, th e accretion rate is reduced to 
values of the order of = .1 x M disc dArtvmowicz & Lubowll 199*61 ; 
iLubow & D'Ange"loll2006h . Here we have considered both the case 
described in Eq. Q, appropriately scaled down for brown dwarfs, 
and the case where accretion is further reduced due to gap open- 
ing, and found that this had little effect on the final distribution of 
planet masses and semi-major axes. Indeed, the details of gas ac- 
cretion do not strongly affect our conclusions, since in most cases, 
as discussed below, we observe only very limited gas accretion in 
brown dwarf protoplanets. 

We again acknowledge that these growth models are probably 
oversimplified. The form for r c given in Eq. l(4) neglects the ef- 
fects of core migration on th e accretion rate of planetesimals onto 
the pr otoplanetary core (e.g. lMcNeil etlTE oO 5t iRice & Armitagel 
2003). The treatment of gas accretion giving rise to the expressions 
for Md, and t k _ h in Eqs. l[5] and (6]l is also highly approximate, 
negle cting the effects of opacity d Adams et ai]|2005l ; Ikom a et al.l 
2000) and all the details evidenced by a numeric al treatment of the 
problem dBrvden et al.ll2000l : |Pollack et alj|l996l) . 



2.3 Luminosity Evolution 

The models of llLlI - llL3l treat the central star as a purely static ob- 
ject, with a fixed luminosity equal to that on main-sequence, 
scaling such that L*/L = (M+/M Q ) 114 . For brown dwarfs, a con- 
stant luminosity model would be completely inappropriate, as no 
constant luminosity is ever reached, and we need to take this into 
account. 

We allow stellar luminosity evolution to take place, account- 
ing for thi s using the pre-main-sequence evolution simulations of 
iTout et all dl999l) . Some sample evolutionary tracks are shown in 
Fig. [2] Compared to a static luminosity model, the net effect of the 
inclusion of luminosity evolution is to raise the luminosity at ear- 
lier times. 

It should be noted that the core accretion model considered 
here, with an initial swarm of planetesimals of mass 10 18 g, may 
have required ~ 10 4 yrs for the pla netesimals to form, although the- 
oretical models vary considerably 1 Weidenschilli ng & Cuzzifl993l : 
ISupulver et al.ll997l : lRice et alj2004ll2006h . Hence we set the zero 
of time in our model to correspond to 10 4 yrs in Fig|2] 



2.4 Planetary Migration 

Protoplanetary cores embedded within a gaseous disc interact with 
the gas and are thus subject to gravitational torques which cause 
them to migrate away from their initial position. For small core 
masses the interaction is linear and the corresponding migration 
regime is called Type I, while for larger core masses (typically 
above 1 Jupiter mass, for planets forming around a solar type star), 
a gap is opene d in the dis c and the resulting migration regime is 
called Type II dWardll 19971) . 

Type I mi gration is important for a num ber of reasons: (i) It 
is well known dPapaloizou & Ter quem 2006) that the inclusion of 
Type-I migration in models of planet formation around solar-mass 
stars completely prevents the formatio n of massive gas giants (un - 
less the disc mass is sufficiently small iThommes & Murrav 2006), 
as all the planetary cores rapidly migrate into the central star, and 



4 M. J. Payne and G. Lodato 




2 4 6 8 10 

log(t/yrs) 

Figure 2. Luminosity evolution profiles from Tout et al. 1 1999) for (top-to- 
bottom) M*/Af Q = 1.5, 1.0, 0.6, 0.4, 0.2, 0.084 and 0.05 respectively. 



th us Type I migration is general ly neglected. However, simulations 
by iNel son & PapaloizQtj J2004T) suggest that, in a turbulent disc, 
Type-I migration may actually follow a random walk, potentially 
providing a mechanism to ameliorate this problem, (ii) It may help 
to explain the presence of intermediate mass planets at small semi- 
major axes, by providing an additional migratory mechanism. In 
this paper, we ha ve run mode ls both including and neglecting the 
effect of standard dWardll 19971) Type-I migration. 

Type-II migration occurs when the planet becomes massive 
enough to open a gap in the disc, the subsequent orbital evolution 
of the pl anet becoming firmly cou pled to the viscous evolution of 
the disc (Lin & Papaloizou 1985). A gap is opened up when the 
planetary mass exceeds 




where we have use d the a-parameterizat i on for the effective viscos- 
ity, v = aH 2 n K , of IShakura & Sunvaevl d 19731) a nd we use c 2 oc T. 

However, hydrod ynamical simulations dLin & Papaloizoul 
ll985l:llda&Lin| [2004a) suggest that the onset of Type-II migra- 
tion occurs somewhat later, when the planet has grown to a larger 
mass. Type-II migration is therefore taken to begin when the planet 
exceeds 10 times the gap opening mass. 

In general, in any viscously evolving gaseous disc, the in- 
ner parts move in and accrete, while the outer parts move out 
and take up the angular mom entum lost by the accretion material 
dLvnden-Bell & P ringle 1974). The direction of Type II migration 
is therefore determined by the position of the planetary core with 
respect to the radius which marks the transition from the inner to 
the outer d isc, which we obtain from th e appropriate self-similar 
solution of lLvnden- Bell&Pringleldl974l) . 



2.5 Methodology 

Given the above semi-analytic prescription for protoplanetary 
growth, if we specify M t , k, n, f , T disc and the initial semi-major 
axis dj, we can then allow the planet to grow from an initial mass of 
10 22 g. If we then perform a Monte-Carlo simulation, allowing the 
parameters f , T disc and a, to vary, we can build up a picture of the 
potential for planet formation for a wide range of protoplanetary 
discs around a wide range of stars. 

Because the initial size distribution is unknown for brown 
dwarf disks, the dissipation timescale is also unknown. At- 
tempts at estimating this timescale have been made by 
I Alexan der & Armitagg d2006l) based on the M - M relation 
dNatta et al]|2004l) . and they argue that t ~ M~ l . However, it is 
not clear whether such observational scalings really reflect some 
intrinsic properties of the initial protostellar dis c population, nor 
how m uch are they affected by selection effects dClarke & Pringld 
2006), so we have kept the same timescale for solar mass stars 
and brown dwarfs (but we d o allow for a r elatively large spread 
in timescale, as sugg ested bv lClarkdd2007l) ). We therefore allow 
Tdisc to be evenly distributed in logarithmic scale over 10 6 - 10 7 yrs. 

We take f to be log-normally distributed (centred on log / = 
and with variance equal to 1) and then demand that the disc mass 
within 100AU be less than 0.35M*, to ensure that the disc is grav- 
itationally stable out to 100AU. We thus end up with a distribution 
of disc masses such that the maximum disc/stellar mass fraction is 
^ 0.35, with an avera ge significantly below this value. Note that 
Rodriguez et al. (2005) have observed protostellar disc masses as 
high as 0.35M*, indicating that the above upper limit is reasonable. 
Finally, we distribute the initial orbital semi-major axes evenly in 
log scale over 0. 1 to 100 AU and terminate all migration at O.O4AL0 



3 RESULTS 

Before discussing our results on planet formation around brown 
dwarfs, we initi ally te st ou r model by reproducing some of the 
main results of IlLlI - llL3l for planets around Solar-type stars. 
Our standard model assumes a disc profile with S oc aT il2 and 
fo oc (M^/Mq) and considers only Type-II migration. We then ex- 
tend our standard model down to the Brown Dwarf regime, adding 
in also the effect of luminosity evolution (section [3751 . Finally, we 
investigate the impact of varying our main parameters on the po- 
tential for planet formation ( section |Tl4t . 

3.1 Planet Formation around Solar-Mass Stars 

Fixing M+ = M B , we reproduce in Fig. |3]the Mass - Semi-Major 
axis distribution as found in IL1, fig. 12c. The model uses our stan- 
dard values such that k = -3/2 in E oc (a/lAU)* and the distribu- 
tion of values of fo is such that the maximum disc mass is 0.35M,, 
with an average of 0.05M*. We use the following nomenclature, 
approximately equivalent to that of llLll to label the planets: 

• Gaseous: M g > WM C ; 

• Intermediate: M c < M g < 10M r ; 

2 While in general, for solar-type stars, this radius corresponds to the stel- 
lar magnetospheric radius, there is no strong reason for assuming that the 
termination radius in brown-dwarf discs would be the same, but as we show 
later, we find little or no migration around brown dwarfs, so its exact value 
is effectively unimportant. 
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Figure 3. Planet masses > iM s for M+ = M Q with static luminosity. 
The different symbols respectively label: Open Circles - "Gaseous" (M g > 
10M C ); Filled Triangles - "Intermediate Composition" (10M f > M g > M c )\ 
Crosses - "Rocky" (M g < M c ). 



• Rocky: M g < M c ; 

where M g is the final gas mass and M c is the core mass, such that 
the total planet mass is M = M g + M c . 

The model distribution presented in Fig. [3] faithfully repro- 
duces that of IlLlL fig. 12c: They both display the characteristic 
"planet desert" for a < 4AU and 10M e < M < lOOMg, and 
both have a maximum gas-giant mass of ~ 3000M ffi . In addition, 
the distribution of planet types within the diagram is very simi- 
lar (although the precise definition of planet-type-labels may differ 
slightly). 

Note the high density of massive planets at a = 0.04AU caused 
by the inward Type-II migration being arbitrarily halted at this dis- 
tance. 

Note also that below the IM S line, there is a vast tail of small- 
mass planets produced by the simulations; planets which grew only 
very slowly and/or were growing in a disc whose mass was very 
small. 



3.2 Planet Formation around Intermediate-Mass Stars 

The results of Fig [3] are in good general agreement with observed 
extra-solar planets, the exception being the problematic observa- 
tions of intermediate mass planets at small semi-major axes, 
note that their results indicate that close-in, intermediate mass plan- 
ets cannot form around F,G and K dwarfs, but may form around 
M dwarfs and note that the formation of dynamically isolated low 
mass planets around such stars would likely be attributable to Type 
I migration or the effect of sweeping secular or mean-motion reso- 
nances. 

In Fig [4] we focus on the observed planetary systems around 
Gl-581 and Gliese 876, both of which have M t ~ 0.3 M Q . We thus 
plot the observed planets on the same scale as our numerical results, 
using models that do not include Type I migration. We consistently 
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Figure 4. Planet masses > lM e for M t = 0.3M o with only Type II migra- 
tion present. Labels are as in Fig. [3] with the addition of observed planets 
plotted as filled squares. 



find it very diffic ult to produce clo se in, intermediate mass planets 
(such as Gl-581-c ludrv et alj2007h . if we do not include some form 
of Type I migration in the process. 

However, as mentioned above, the simple introduction of 
the fast Type I migration (Ward 1997) would generally pre- 
clude the formation of the high mass (~ 100M e ) planets ob- 
served in these systems, unless the m igration rate is significantly 
slower than the standard an alytic form (Nelson & Papaloizo ul2004l ; 
iThommes & Murray||2006f) . It thus appear that, while some addi- 
tional form of planetary core migration is needed in order to repro- 
duce the above-mentioned observations, its rate should be substan- 
tially reduced with respect to the simple analytical estimates. In the 
rest of the paper we will therefore consider the two extreme cases 
in which we neglect completely Type I migration and in which we 
include standard Type I migration. 



3.3 Planet Formation around Brown Dwarfs 

If we now allow the central stellar mass M+ to vary down below 
0.084M o , we can to investigate the potential for planet formation 
in the brown dwarf regime. As noted in section [2731 this extension 
requires the implementation of a luminosity evolution model for 
brown dwarfs. 

In Fig.|5]we plot the mass - semi-major axis distribution result- 
ing from our numerical simulation of the planet formation process 
around brown dwarfs of mass 0.05M o - note the change in scale 
compared to Fig. [3] We keep the standard disc profile of section |3~T| 
and assume that the disc surface density scales as Mj (i.e., we as- 
sume n = 2 in Equation $2%). This effectively makes the 0.35M* 
cut-off unnecessary, as the maximum disc mass is now 0.21M* 
(which corresponds to 10.5Mj upiCer ), while the average disc mass 
is0.003Af*(0.15M Jupiter ). 

Clearly the outcome of the planet formation process around 
brown dwarfs is radically different to that for G-dwarfs. Of partic- 
ular importance is the fact that all of the planets produced are of 
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low mass - we do not see any significant number of planets formed 
above IM S , the most massive ones corresponding to the high tail of 
the disc mass distribution - and are of an icy constitution. In none 
of the simulation runs did any planet grow massive enough for any 
significant gas accretion take place. This means that the distribution 
of planetary masses is effectively limited by the surface density of 
the solid disc. The radial extent of the planet population is at least a 
factor of ten smaller than around solar-mass stars and we find that 
the peak of the distribution shifts to significantly smaller radii as 
the stellar mass decreases, reflecting the fact that disc densities at 
any given radii are lower around brown dwarfs. 

Whilst the lack of gas accretion limits the predicted mass of 
the planets, it also has implications for planetary migration. If a 
planet does not accrete gas and remains relatively low mass, it is 
highly unlikely to open up a gap in the gaseous disc, and hence will 
not under-go Type-II migration. This means that the only migra- 
tion mechanism relevant for the brown dwarf regime will be Type-I 
(which has not been included here). Thus, any observations of plan- 
ets around brown dwarfs could shed light directly on the effects of 
Type-I migration and hence illuminate its likely character in planet 
formation around all higher mass stars. 

Note also the paucity of planets plotted in Fig. [5] indicating 
that Earth-like planet formation is highly unlikely in this case. In- 
deed, we found a planet with mass M > 0.3M e in only 0.035% of 
the simulations that we have performed. 

3.4 Impact of Parameter Variation on Planet Formation 

3.4.1 Surface Density Dependence on Semi-Major Axis 

It was assumed in section [331 that the disc profile varied as Z oc 
a~ 3/2 . In this section we explore the effect of having a shallower 
density profile, with E oc a -1 , while keeping the scaling with stellar 
mass as (M^/Mq) 2 , and allowing only Type II migration. We find 
that at both solar mass and lower masses of M+ = 0.05M o , the 
effect of moving to k = -1 is significant. 



At = M , moving to &k = -1 model leaves the average 
and maximum disc masses unchanged from the standard model in 
section [37X1 However, a comparison between Fig. [6] (which shows 
the planet mass - semi-major axis distribution in this case) and Fig. 
[3] shows that for semi-major axes < 10AU, the change in model 
gives very similar results for the distribution of Jupiter-mass plan- 
ets, but suppresses the production of terrestrial planets. Fig.[7]shows 
the average and maximum planetary masses at different radii for 
the various cases. This shows that both the maximum and average 
planet mass inside 10AU (which are dominated by the gas giants) 
are similar for these models, but beyond 10AU, the effect of the 
k = - 1 model is to increase both the average and maximum plant 
masses. This result can be simply understood based on the fact that, 
for a given disc mass, with a shallower surface density profile the 
disc dens ity is decreased at sma ll radii and increased at large radii. 
Note fhat lRavmond et al] d2005t) looked at terrestrial planet forma- 
tion at small radii and found it was suppressed in discs with shal- 
lower profiles, in agreement with our result. 

Let us now turn to analyse the case of brown dwarf. Fig. [8] 
shows the planet mass - semi-major axis distribution for k = - 1 
(cf. Fig. [5j- The effect of moving to a J: = -1 model, similarly 
to the case of higher stellar mass, is to suppress the formation of 
rocky planets. This is further evidenced in Fig. [9] which shows a 
comparison of the maximum and average planet mass at different 
radii for the two cases k = -1 (solid line) and k = -3/2 (dashed 
line). Unlike the case of Solar-type stars, these histograms are now 
dominated by rocky planets, and indeed moving to a shallower sur- 
face density profiles results in a reduction of both the maximum 
and average mass planet formed, preventing the formation of any 
planets with masses greater than lM e . 

3.4.2 Surface Density Dependence on Stellar Mass 

Disc masses estimates for brown dwarfs are very un certain and only 
known for a handful of objects ( Scholz et al. 2006j). In our standard 
model, (section [3~3t we have assumed that protostellar disc masses 



Planet Formation around Brown Dwarfs 7 



2 - 



CUD 
O 




Figure 7. Changing the slope of the surface density profile at = 1M Q . 
The thin lines indicate average masses, the thick maximum masses, whilst 
solid lines indicate k = — 1, dotted lines k = -3/2. 
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Figure 9. Changing the slope of the surface density profile at M t = 
0.05M G . The thin lines indicate average masses, the thick maximum 
masses, whilst solid lines indicate k = -1, dotted lines k = -3/2. 
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Figure 8. Planet Masses versus semi-major axis for M+ = 0.05M Q with 
k = — 1. Labels are as in Fig. [3] 



scale with stellar mass as Ml , leading to typical disc masses of a 
tenth of a Jupiter mass in the brown dwarf regime. Here, we exam- 
ine the effect of allowing for relatively more massive discs around 
brown dwarfs, and we thus consider the case where the disc surface 
density scales linearly with stellar mass. 

We plot in Fig.[l0]the resulting planet mass versus semi-major 
axis distribution in this case, for the two cases where k = -3/2 
(steeper density profile; left panel) and k = -1 (shallower density 
profile; right panel) and collate in Fig.Qjjthe histograms of average 



and maximum planet masses formed around brown dwarfs for the 
various models. 

In this case, typical disc masses are of the order of a few 
Jupiter masses and the enhanced disc mass is reflected in a larger 
fraction of objects being able to form Earth-like planets. Indeed, 
planets with mass larger than 0.3M e are formed now in «10% of 
the cases for k = -3/2 and in re 2% of the cases for k = -1. 

In this case, the cores are growing faster, but still none of 
them grow quickly enough to start accreting gas before it dissi- 
pates. Thus, despite the increased potential for planet formation, 
no object is able to accrete a significant gaseous envelope, and the 
maximum planet mass is only slightly enhanced, the maximum now 
being 5M m , in the most favourable case (that is, k = -3/2). 

The maximum planetary mass obtained is also slightly depen- 
dent on the upper limit on the disc mass in order for the disc to be 
gravitationally stable. For example, if we only allowed discs less 
massive than 10% of the central brown dwarf mass, the maximum 
resultant planetary mass would be decreased to 1 M e and it would 
be raised to 9M S , with some very limited gas accretion occurring, 
if the discs are allowed to be half as massive as the brown dwarf. 

Clearly, an important parameter in the determination of the fi- 
nal planet population is the disc surface density. For a given total 
disc mass, a more compact disc, with an outer radius significantly 
smaller than the 100 AU assumed above, would have a larger sur- 
face density and therefore an enhanced chance of forming planets. 
We have thus re-run the most favourable case for planet formation 
(that is, the one with n = 1 and k = -3/2), assuming the same dis- 
tribution of total disc mass, but with significantly smaller disc sizes 
(this, in turn, implies increasing the factor fo in Eqs. (1) and (2)). 
The results are shown in Fig. [12] for an outer disc radius R out = 50 
AU (left panel) and 20 AU (right panel). While the reduction of 
/? ut to 50 AU only results in a slightly increased maximum planet 
mass, in the case of R out = 20 AU we do see the appearance of a 
number of gas giants in a limited number of cases. However, we 
regard these cases as extremely unlikely, since they correspond to 
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Figure 10. Planet Masses versus semi-major axis for Brown Dwarf with M+ = 0.05M o , n = 1 and k = -3/2 (left panel) and k = -1 (right panel). Labels are 
as in Fig. [3] 
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Figure 11. Results for planet formation around brown dwarfs. The left panel shows the average planet masses while the right panel shows the maximum planet 
mass in the various models considered. Solid Line - n = 1 and k = -3/2; Dotted Line - n = 1 and k = — 1; Dashed Line - n = 2 and k = -3/2; Dot/Dash Line - 
n = 2 and k = — 1 . 



the upper tail of the disc mass distribution, with disc masses of the 
order of several Jupiter masses, confined within a region as small 
as 20 AU. 



3.4.3 Introduction of luminosity evolution 

As mentioned above, the models that we use for planet formation 
around brown dwarfs also include the effect of the evolution of the 
(sub)-stellar luminosity with time. The stellar luminosity has two 
main effects on our models: (a) it determines the temperature of 
the gas and (b) it sets the location of the snow-line, inside which 
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Figure 12. Planet mass versus semi-major axis for M 4 = 0.05M Q , n = 1 and k = —3/2, with outer disc radius equal to 50 AU (left panel) and 20 AU (right 
panel). While in the case of R out = 50 AU, the final distribution of planet mass is very similar to the case of R out = 100 AU, some gas giants appear in the 
case R out = 20 AU. These latter cases, however, correspond to the very unlikely case of a massive disc (M L \[ SC > a few Jupiter masses) with a very small outer 
radius. 



the solid density is significantly reduced (we neglect here any other 
possible effect related to changes in the efficiency of angular mo- 
mentum transport within the discs). Since there is little or no gas 
accretion on proto-planets around brown dwarfs in our model, and 
no planet becomes massive enough to open up a gap and start mi- 
grating, the first of the above effects actually has essentially no im- 
pact on the final outcome. The second effect is indeed present, but 
it turns out not to affect deeply the final results, and simulations 
performed without taking into account the luminosity evolution of 
the brown dwarfs only show relatively small differences. 



in this case is relatively longer. For the n = 1 k = -3/2 model (see 
Fig. 1 131>, the maximum planet mass at ~ 1AU is very slightly re- 
duced, whilst the average value at this radius (and the efficiency of 
the planet formation process) is reduced as a large number of the 
cores are swept into the inner boundary. 

Since Type-II migration effectively never occurs around 
brown dwarfs, any observational evidence of a secondary 
peak/excess in the radial distribution of planetary semi-major axes 
inside the main distribution peak at ~ 1AU would provide a clean, 
direct insight into the nature and extent of Type-I migration. 



3.4.4 Introduction of planetary migration 



As IIL3I have already noted, since low-mass cores do not undergo 
Type II migration, any detection of close-in super-Earths around 
low mass sta rs, as in the case o f Gl-581-c (a 5.5M e planet orbiting 
a 03M e star. lUdrv et alj|2007l) . would need the inclusion of other 
types if migration. In our models of planetary formation around 
brown dwarfs in which Type I migration is neglected, we always 
predict a peak in the distribution at 1AU, irrespective of the other 
parameter variations in the model. 

As discussed in section [3~2l we now consider the addition of 
a form of Type I migration to our model. In p articular, we consider 
the model for Type-I migration as derived by iTanaka et al.l ( 2002) 
for a 3-dimensional isothermal gaseous disc, where the migration 
rate is of the form 



Tl 



1 x 10 5 



3/2 



2 x 10 3 g cnr 



yrs, 



(9) 



We find that in general the effects of Type I migration around brown 
dwarfs are less significant than around solar type stars. This is 
mostly due to the fact that planetary cores around brown dwarfs 
grow at a much slower rate, and therefore the migration timescale 



4 DISCUSSION AND CONCLUSIONS 

We have extended standard core-accretion models for planet for- 
mation down to the brown dwarf regime and investigated the po- 
tential for planet formation in protoplanetary discs around objects 
with masses ~ 0.05M G . We looked at the impact on the predicted 
planetary mass - semi-major distribution of varying fundamental 
parameters, like the total disc mass and i ts radial distribution. 

In line with previous results (IL1), at higher stellar masses 
(M» ~ 1M , Fig. [3}, we found that the full range of extant planet 
types is produced, including rocky and icy cores up to 10 2 M e , rock 
and ice giants spanning a range 10 1 - 10 3 A/ ffi and gas giants extend- 
ing from 10 2 M ffi up to 10 4 M ffi . In addition, we have found that a 
shallower surface density profile results in a reduced efficiency of 
the formation of rocky planets at radii < 10 AU. 

Our main result is the determination of the likelihood of planet 
formation around brown dwarfs, with mass ~ 0.05 M Q . Our main 
findings can be summarized as follows, (i) Giant planet forma- 
tion is completely inhibited in this case. In none of our standard 
simulations did any planet accrete a significant gaseous envelope. 
This occurs despite the fact that in a few cases relatively massive 
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planets inside this radius would give clear evidence of the nature of 
Type-I migration. 
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Figure 13. Effect of analytic Type I migration on brown dwarfs with M 4 
0.05M o and n = 1, k = -3/2. Labels are as in Fig.fJ] 



cores do form (in principle massive enough to start runaway gas 
accretion), because such massive cores, when they form, only do 
so at a late time, when most of the gaseous component of the disc 
has dissipated. This means that Jupiter mass companions to brown 
dwarfs, as sometimes observed, have to be formed through a dif- 
ferent - binary like - mechanism, that is either core or disc frag- 
mentation. While this had already been discus sed for very wide ex - 
tremely low mass binaries, like 2MASS1207 dLodato et al]|2005h . 
we can now extend this conclusion to closer binaries, (it) The for- 
mation of Earth-like planets is possible even around brown dwarfs, 
and the maximum planet mass that we have found is « 5M ffi . How- 
ever the likelihood of this process is critically dependent on what 
the typical disc masses are in this regime. In particular, if typical 
brown dwarf disc masses are of the order of a few Jupiter masses 
(which is consistent with a linear scaling of disc mass with stel- 
lar mass), then up to 10% of brown dwarfs might possess planets 
with masses > 0.3M ffi , while on the other hand if typical brown 
dwarf disc masses are only a fraction of a Jupiter mass (con- 
sistent with a quadratic scaling of disc mass with stellar mass), 
then even the formation of Earth-like planets would be essentially 
prohibited. We note that current estimates of brown dwarfs disc 
masses are only available for a small number of objects and are 
very close to our boundary between planet-forming and non-planet- 
forming systems, w ith upper limits at the level of a few Jupiter mass 
(Scholz etal.ll2006h . In order to answer the question of the occur- 
rence of planets around brown dwarfs it is therefore essential to 
improve the statistics and the accuracy of disc mass measurements 
in the very low mass regime. 

Finally, we note that the lack of any gas accretion in the planet 
formation process around brown dwarfs and the concomitant lack 
of any Type-II migration means that the radial distribution of the 
planets will be determined purely by the initial mass distribution 
in the disc and the subsequent effects of Type-I migration. Given 
that the underlying disc models tend to produce the most massive 
planets at a ~ 1AU, irrespective of the specific model, then any 
excess/secondary peak (or lack thereof) in the number of detected 
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